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r-j ■ We investigate possibilities for detecting heavy Majorana neutrinos (iV's) in 

Q^h, e + e~ at LEP200 and future Linear Colliders. We concentrate on the processes 

CX| where the pairs of intermediate heavy iV's produce a clear signal of total lepton 

number violation (e + e~ — * NN — * W + l~W + l'~). Such a signal is not possible 
if the heavy neutrinos are of Dirac nature. Our approach is general in the 
sense that the intermediate iV's can be either on-shell or off-shell. Discussion 
of the relative numerical importance of the s and the t+u channels of the NN 
production is also included. 
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There has been a significant amount of activity in the high energy physics community 
towards discerning the nature of the neutrino sector. A basic question is: Are neutrinos Dirac 
or Majorana particles? If there are no right-handed currents, then it is virtually impossible to 
discern the nature of the light neutrinos Jl|. If there are heavy neutrinos (M~ 10 2 GeV), then 
present and future experiments offer a realistic prospect of establishing their nature. The 
production cross section of heavy Majorana neutrinos (iV's), mostly via the e + e~ collisions, 
has therefore been investigated in the past 0. Most of these works have been done within 
specific (classes of) models; and it has been assumed that the center-of-mass (CMS) energy 
y/s in the process is high enough for the production of on-shell (OS) heavy iV's. The effects of 
the off-shell (nOS) iV's have been ignored. Moreover, to our knowledge, various distributions 
of iV's decay products [N — > W ± P(—> jets + £ T )], which are produced in the full reaction 
and which can actually be detected, have not been investigated in a quantitative way. The 
main reason for this was that the expressions for invariant amplitudes with two on-shell iV's 
apparently do not allow a straightforward calculation of such distributions. We note that 
the detection of events for the reactions e + e~ — ► NN — ► W £^W i'^ 1 — > jets + £ T £' T , which 
violate the total lepton number, would be a clear signal of the Majorana character of the 
intermediate neutrinos. 

We present some results of calculations for the afore-mentioned reactions. We do not 
restrict ourselves to any specific (classes of) models. In contrast to the available literature, 
our approach allows us to account also for the effects of off-shell (nOS) intermediate iV's on 
the cross sections a. This enables us to investigate deviations from the previously known cr's, 
in the "20S" kinematic region (y/s>2M >2M\y) where both intermediate iV's can, but need 
not, be on-shell - these deviations are termed "finite width effects." Our approach allows us 
to calculate the cr's, and various distributions, also in the "10S" region (2M> y/s> M+Mw) 
where at most one intermediate N can (but need not) be on-shell, and in the "nOS" region 
(M+Mjy > y/s) where both iV's always have to be off-shell. Our approach makes possible a 
straightforward calculation of various distributions of final particles. As an illustration, we 
include an angular distribution of the final state leptons £, £'. 

We start with rather general Lagrangian densities for the couplings of N with Z, and for 
the coupling of iV with W and light leptons £j (£j~ = e~, £^~=/i~, £ 3 = t~): 

Cnnz(x) = - 9 — A NZ N(x)^ 5 N(x)Z fl , (1) 

4 COS t>w 
3 n(j) _ 

C mw (x) ^-J2^£ j (x)ri-N(x)W-+h.c, (2) 
i=i 2 ^ 2 

where 7_ = (1 — 75); B^'s are, at first, free parameters; g and 9w are the standard SU(2)l 
gauge coupling parameter and the Weinberg angle, respectively. The vector part is absent 
in (TJ) because iV's are Majorana. The right-handed parts were neglected in (fj). The other 
relevant coupling is e + e~Z which we consider to be the one of the Standard Model (SM). We 
also set A NZ = —^, i-e., by replacing, in the SM density for uuZ, the massless Dirac neutrino 
v by the (heavy) Majorana neutrino N. These choices would suggest that the considered iV 
is made up primarily of a sequential neutrino with the standard SU(2)lxII(1)y assignments. 
However, these choices may also represent an approximation to other scenarios (cf. @, @, 
0). Further, parameters in (||]) will affect the final results only via the combinations 
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m=i4 1} r, H=j:\B<!>\*. (3) 

i=i 

We restrict ourselves to the afore-mentioned reactions e + e _ — > NN — > W^ifW^iJl—* 
jets + £^£j) with light leptons (£1 = e, £2 = /•*, ^3 = t). They involve the s and the t + u 
(shortly: tu) channel - cf. Fig. [I]. 

For the calculation of the invariant amplitude Aifi (shortly: M. P|) for various channels, 
we used the 4-component spinors u^ a \q) = u(qa) and t> ^(5) = f (qa) as defined in 0, but 
with the normalization convention as given in ||. For the quark spin components we use 
the notation: 5 = 1, 2 a — 2, 1. In the s-channel, the resulting amplitude is 

[ Pjy (p#iu ) -P/v (P£p™ )u{p t a t )C^ {pwPwPiPi ) ( 1 + 75 ) v {pifin ) + {p w ^Pw)} ■ (4 ) 

Here, we use notations of Fig. [I]; s— (p+p) 2 ; Mz and Tz are the mass and the total decay 
width of Z; A { + ] and are the vector and axial- vector coupling parameters of the e + e Z 
coupling of the SM, respectively (A+ ^ =4 sin 2 Qw~ 1, A_ =— 1). C^'s are: 

&(p w p wP m) = t [&+jOt* + -ftft+tj] t , (5) 

where $ = 9^7^; = s^(p w ) and e„ = £^(p M ) are the real polarization vectors J4| of the 
final W's, with polarizations A, A = l, 2, 3. M is the mass of iV's; is: 

A« = g 2 B^B^A NZ i\ M /(128 cos 2 , (6) 

where Am is the phase factor in the Fourier decomposition of the Majorana field N(x) (cf. ||; 
|Am| 2 = 1)- Pn in @ is the (scalar) denominator of the propagator of iV 

P N (PiP w ) = l/[( P i+p w ) 2 -M 2 + iMT N ] , (7) 

where T N is the total decay width of N. 

All the combinations of attaching the four final particles to the two iV's are accounted for 
in amplitude (|j). Further, (Q) had originally four instead of two terms in the curly brackets; 
however, two of them were reduced to the other two, by using the general identities 

-rfu(qaY = (-l) a v(qa), -ij 2 v(qa)* = (-lfu(qa), (8) 

and (4P) T = ~ 7°7 2 W)7 2 7°> where the Dirac basis and the conventions of || are used for 
7 A "s. Using (g), we can further rewrite (||) into a form involving v(p£ai) and u(ptat) instead 
of u(peae) and v{p[ai). The complex conjugate (c.c.) of this alternative form and of the 
form (|]) are needed to calculate later the s-tu interference term of the full (|.M| 2 ), where 
(...) stands for summation over the final and average over the initial polarizations. For the 
s-s part of (|.M| 2 ), the form (f|) is used. 

The tw-channel amplitude turns out to be 
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iM (tu) = ±MAWp N (p ePw )P N (p e p w ) 5 
~ [(p-pe- Pw ) 2 -M^ + iT w M w } [ ' 

| u(p e a e )A(p w pe)Yl~u(pa) v(pa)^ u ^ + v(pia t ) + 

+ [u(Pi a t)fa- u (P a )} \v(pa)A(p w pe)j+v(p e ae)\ } + ..., (9) 

J-VLyy 

where A(p w pi) = (fy w + 2p e ■ e), A(p w pe) = (p w $+2p l ■ e), j ± = (1±7 5 ). The dots at the 
end of (|9|) stand for three analogous terms, obtained from the above explicit expression by 
replacements: I. (p w , e) <-> (p w , e); II. (pi, ai) <-»• (p^, o^) and overall factor (—1); III. combined 
replacements I. and II. in @ is 

A^=g 4 \B^B^B^i\ M /6A. (10) 

As in the s-channel case, we can reexpress any of the terms in MS tw > in alternative forms, 
by applying transformations @ - e.g., if we want to use, in scalar expressions in square 
brackets of @, u(pa) and u(picei) instead of v(pa) and v(pia). Such transformations are 
convenient when we calculate (|.M| 2 ) = (\A4^+Ai^\ 2 ). Then we can always end up with 
traces involving u(q, fl)u(q, (3) and/or v(q, j3)v(q, 0)=4 (q = P,P, ■ ■ ■] /3 = a, a, . . . = 1, 2). 

The integrand (|.M| 2 ) is long - the s-tu and (even more so) tu-tu parts extend over tens 
of pages when printed out. Numerical integration of (|.M| 2 ) over (parts of) the final phase 
space leads to the cross sections. This general (nOS) program, as mentioned, accounts for 
the effects of off-shell and on-shell iV's. 

The input were values of y/s, M, HI and H [cf. Eqs. (0)-©]- HI measures the eWN 
coupling and affects the tw-amplitude (oc HI). H affects the total (|.M| 2 ) which is then 
formally ocH 2 (if HI is fixed). In (|.M| 2 ), we average over the initial (a, a) , and sum over 
the final polarizations (A, A; ae,ae) and over the flavors (i, j = 1, 2, 3) of the two final light 
leptons. In the general (nOS) expression, an additional factor 1/4 is included in (|.M| 2 ) to 
avoid double-counting of the two iy + 's and of the final leptons when integrating over the 
phase space. 

By the same described methods we also calculated M. and (|.M| 2 ) when one N, or both 
N's, are explicitly put on-shell (10S, 20S expressions, respectively). The 10S (|A^| 2 ), for 
the sum of reactions e + e~ — > NN Q s — » W + £^Nos (i— 1,2,3), was then multiplied by the 
branching ratio Br for the sum of the decay modes iVos — ► (j = 1,2,3); the 20S 

(\M\ 2 ) for e + e~ — > N 0S N s was multiplied by (Br) 2 . 

Tn, appearing in (|^), was determined at the tree level, assuming that the only (dominant) 
decay modes are N^W ± £J (J = 1,2,3) [=>T N ocH]. Then Br = 1/2. 

Numerical calculations were performed in various kinematic regions (nOS, 10S, 20S 
regions) with the general (nOS) expression [cf. (171)]. In the 10S and 20S regions, the 10S 
expression was also used. In the 20S region, the 20S expression was also used. Results are 
depicted in Figs. ||-|| The IV-parameter H was set H = l in all these Figures. 

Figure |2| shows the M-dependence of the cross section a, at fixed t/s. The difference 
between the results of the general (nOS) and the 20S program, for = 300 GeV, is less 
then ten percent over most of the 20S kinematic region (My/ < M < y/s/2), except near 
the threshold M^^fsj2 where the results of the nOS program are significantly higher. The 
difference between the results of the 10S and 20S programs is less than three percent in 
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most of the 20S kinematic region. However, in the 10S region (y/s/2<M < y/s— Myy), the 
results of the 10S program are usually by several factors lower than those of the general 
nOS program, except near the threshold MfHy/s/2 where they differ only little. All these 
differences are in general less pronounced when the tu-channel is excluded (if 1 = 0.0). For the 
chosen tii-strength value HI = 0.25, the contributions of the tu-tu channel are at least by one 
order of magnitude larger than those of the s-s channel. Each curve has two slope increases: 
at M^iy/s/2 and at M ~ y/s—Mw (onset of the 20S, 10S kinematic region, respectively). If 
we take the integrated luminosity at LEP200 (^=200 GeV) to be 500 pb -1 , Fig. || implies 
that the maximal number of events cannot exceed 17 and 112 if HI — 0.0, 0.25, respectively 
and we assume M > 85 GeV. 

In Fig. ^| we show the v^s-dependence of cr, at fixed M. Most of the remarks about Fig. |2| 
apply also to Fig. |^. The differences between the results of the nOS and 20S programs 
("finite width effects") are very significant when M = 200 GeV (and HI = 0.25, # = 1.0), 
because Pjv = 4.8 GeV is large then (r7v/M^2.5%; for M = 150 GeV: T N /M ^1.1%). 

Our general nOS program can be applied also to calculation of various distributions in the 
process. As an illustration, we present in Fig. [| an angular distribution of the final leptons 
£~, t'~ . The corresponding total cross sections are a = 0.280 pb, 0.007 pb, for M — 200, 255 
GeV, and the kinematic regions are 20S, 10S, respectively. If linear colliders at ^^=500 
GeV achieve the integrated luminosity of 10 4 pb -1 , and most of the final states W + W + £~£~ 
can be identified, then these cr's will correspond to 2800 and 70 events, respectively. 

In the 20S kinematic region, the cr's and distributions as given numerically by the general 
(nOS) program depend on H weakly. Parameter H (note: TnocH) is responsible in the 20S 
region for the deviation of the full a from the pure 20S cr. In the 10S region, ii-dependence 
of the full cr becomes quite strong (approximately croc if), and in the nOS region even more 
so (ocH 2 ). In Figs. we chose H = l. 

In large classes of models, in which heavy neutrinos are sequential or have exotic SU(2)x 
U(l) assignments, HI and H of (|3|) are severely restricted by available experimental data 
(LEP and low-energy data) ||[7]]: HI < 0.016, H < 0.122. In principle, in certain models these 
restrictions can be avoided. Taking into account the restriction iiT < 0.016, the influence 
of the tw-channel is much weaker than for HI = 0.25, but the s-tu interference term still 
increases the cross section significantly above the pure s-s contributions (cf. Table |). The 
restriction H < 0.122 would imply that the displayed off-shell ("finite width") effects in the 
20S region would be weaker than in the H=l case, and that the displayed cr's in the 10S 
and nOS regions would be reduced, approximately by factors H and H 2 , respectively. The 
numbers in Table p] were obtained from the general (nOS) program, except in the case of 
M = 85-95 GeV when the 20S program was used. 

To summarize, we calculated cross sections for e + e~ —>■ NN — > W + £~W + £'~ where iV's 
are Majorana neutrinos (mass M ~ 10 2 GeV) and £, £' are light leptons e, fi, r. In contrast 
to the calculations available so far, we included the effects of the off-shell intermediate iV's. 
These effects were significant even when y/s > 2M (> 2Mw)- They are more pronounced 
when the t + u (shortly: tu) channel contributions are significant. The number of reaction 
events at LEP200 (y/s = 200 GeV; integrated luminosity ^500 pb -1 ) would be low (<20, 
for M > 85 GeV) if the strength of the tw-channel is restricted by available experiments 
and by confining ourselves to certain classes of models where the heavy neutrinos are either 
sequential or have exotic SU(2) x U(l) assignments (i.e., when iil < 0.016). Numbers of 
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such events in general significantly increase at linear colliders = 500 GeV; integrated 
luminosity 10 4 pb _1 ), and may be significant even when 2M > y/s. Further, our approach 
allows us to calculate numerically various distributions of the final particles. 

We ignored the questions connected with the experimental difficulties of detecting the 
discussed process unambiguously. In particular, there are problems connected with identifi- 
cation of the (on-shell) Ws and r's. Further, we ignored the possibility M < Myy (GC thanks 
Amitava Datta on that point) - however, additional problems arise in the identification of 
the process since the two W + 's are then intermediate off-shell. 

GC thanks A. Datta for helpful discussions and D. Schildknecht for financial support 
via the BMBF project no. 332-4005-05-7BI92P toward the end of the work. CSK thanks 
W. Buchmuller for helpful discussions and acknowledges the financial support of Korean 
Research Foundation made in the program year 1997. 
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\Js [GeVJ 


M [GeVJ 


a [pb] 


increase 


500 


200 


0.52 • 10~ 2 (0.77 • 10 -2 ) 


49% 


500 


onn 

200 


fn 1 n — 2 / n '7rr i n — 2\1 

[0.53 • 10 (0.75 • 10 )\ 


[43%J 


500 


255 


0.86 • 10 (1.41 • 10 j 


657o 


500 


255 


[1.14 -10" 5 (1.88 -10" 5 )] 


[65%] 


300 


145 


1.33 -10- 3 (1.69-10- 3 ) 


27% 


300 


155 


0.41 • 10" 4 (0.52 • 10" 4 ) 


27% 


200 


85 


0.34 -10- 1 (0.40 -10- 1 ) 


17% 


200 


95 


0.72 • 10~ 2 (0.84 • 10^3) 


18% 


200 


105 


1.70 -1Q- 5 (1.96 • 10~ 5 ) 


15% 



TABLE I. Values of cross sections a, for various values of y/s and M, and for the iw-strength parameter 
HI = 0, and in (. . .) for HI = 0.016. Given are also the relative increases of a when HI =0^ 0.016. The 
TV-decay parameter H is taken H=l; numbers in [. . .] are for if = 0.122. 
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3+e" -> NN -> W+W+^T = e , ^ 2 = /z, £ 3 = t), as function of 
neutrino mass M, for = 200, 300 GeV and the iu-strengths HI = 0.0, 0.25. Full lines are results of the 
general (nOS) program. Triangles and crosses are results of the 10S and 20S programs, respectively. 
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<s [GeV] 

FIG. 3. Sum of the cross sections for the mentioned reactions, as function of the CMS energy y/s, at 
fixed M = 150, 200 GeV. Again Hl = 0.0, 0.25. Results of various programs (nOS, 10S, 20S) are displayed 
in the kinematic regions where they are applicable. 
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FIG. 4. da/dcos9, where 9 is the CMS angle between £ i; £j. Input values (H, HI, ^fs, M) are denoted. 
The points are results of calculations, and the curves are parabolas (in cos 9) fitted to the points with equal 
weights. Fluctuations are due to the limited statistics of the Monte Carlo integration. 
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